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c/3 I We study the properties of diffusion controlled reactions with traps or reaction sites 

attached to polymer chains or to a star polymer. Using a field theoretical renormalization 
group approach we obtain the scaling of the moments of concentration of diffusing particles 
. near the core of a star polymer and calculate numerical values of the exponents governing 

r^ ! the scaling of the reaction rate. 

O, 
, O, . Chemical reactions between diffusing reactants play an important role. Examples of 

such reactions can be found in different systems, ranging from biological systems to nuclear 

►^ ^ reactors (see e.g. |]I[ and references therein). One more place where these reactions appear 

[*»-. ■ as a limiting stage are aggregation models [0. Of particular interest are reactions between 

^ . reactants of different nature: particles A which freely diffuse in a solution and particles B 

Q ! which are attached to polymer chains immersed in the same solvent P] (the concentration 

of polymers being low enough to allow to neglect the inter-chain interaction). Such a 

Q>^ I process may be considered also as a trapping reaction of n particles of A type and traps 

B: 

A" + 5^0. (1) 
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The reaction rate fc„ of (|I]) in the vicinity of a certain trap on the polymer of size R is 
^ I proportional to the averaged moments of the concentration p of diffusing particles near 

^ ■ this trap. It scales with R as 



X 



K ~< p" >~ {R/l)-'\ (2) 

with / being a characteristic length scale. 

The goal of this paper is to investigate the scaling properties (0) in the vicinity of an 
absorber in form of a star-shaped polymer. We are interested in the scaling properties 
of (0) near the trap placed at the center of a polymer star, provided all polymer chains 
are decorated by traps and particles meeting them are absorbed. As a model we consider 
the polymer chains to be either random walks (RWs) or self -avoiding walks (SAWs). To 
this end we generalize the idea of Cates and Witten |^ to describe both the polymer 
absorber and the random walks of diffusing particles in the same way. In a steady state 
limit the diffusion equation describing particles A is reduced to the Laplace equation. 
In terms of the path integral solution of the Laplace equation one finds that p{r) at 



point r near the absorber is proportional to the number of RWs that end at point r and 
avoid the absorber [§, |^ . The nth power of p(r) is proportional to the nth power of the 
above mentioned number, i.e. it is defined by a partition function of a star of n walks 
[0]. Furthermore, introducing the mutual avoidance conditions between the 'n-walk star' 
and the m-polymer star (representing the absorber) one has to calculate the partition 
function of a co-polymer star consisting of chains of two different species that avoid each 
other. These correspond to the trajectories of diffusing particles (being RW) and the 
absorbing polymer which for the purpose of present study is chosen as a RW or a SAW. 
Making use of the previously developed ||^ theory of copolymer stars and networks and 
mapping the model of co-polymer stars to the appropriate field theory p, ^] one may 
relate the spectrum of exponents (^ to the exponents that define the scaling properties 
of co-polymer stars. In particular, for a co-polymer star consisting of m chains of one 
species and n chains of the other species the number of configurations Z^ scales as [^ : 

Here, the rjmn represent the co-polymer star exponents. The latter have been calculated 
using a field theoretical renormalization group approach |^ and are known in the third 
order of perturbation theory p. By means of a short-chain expansion the set of 
exponents rjmn ® can be related to the exponents Xmn ®, that govern the scaling of the 
rate of reaction (|l]) in the case of an absorber in the form of an m-arm star. Considering 
the absorber to be either a RW star or a SAW star we define the exponents: 

XRW ^ G u) 

\SAW ^ U , U 
mn Imn ' ImQi 



here, we keep notations of [Q for the expressions for exponents r/ defined in different fixed 
points of the renormalization group transformation. Based on these expressions we get 

Table 1: Values of the exponent \^ at c/ = 3 obtained by e-expansion (e) and by fixed 
dimension technique (Sc?). 

~i^ 1 2 3 4 5 6 

n 6 3d e 3d 6 3d e 3d e 3d e 3d 



1 


0.56 0.58 1.00 1.00 1.33 


1.35 


1.63 


1.69 


1.88 


1.98 


2.10 


2.24 


2 


1.77 1.81 2.45 


2.53 


3.01 


3.17 


3.51 


3.75 


3.95 


4.28 


3 


3.38 


3.57 


4.21 


4.50 


4.94 


5.36 


5.62 


6.15 


4 






5.27 


5.71 


6.24 


6.84 


7.12 


7.90 


5 










7.42 


8.24 


8.50 


9.54 


6 














9.78 


11.07 



the result for the exponents in the form of an e-expansion series: 

e e'^ e 

X^ie) = mn- + mn{3 — m — n) — [1 — (m -|- n -|- 3^(3) — 3)-], (5) 

2 8 2 



e e^ 

Xi^^ie) = 3mn--mn(18n + A2m-91) + mra [492 m^ - 2290m + 2463 + 

8 256 

108 n^ - 1050 n + 540 nC(3) + 1188 mC(3) + 504 mn - 2652 C{3)]——.{6) 

4096 

Here, e = 4 — c?, d is the space dimension, and ({3) ~ 1.202 is the Riemann zeta function. 
Similarly, we obtain the expansions for (^) in form of a series at fixed d = 3 dimensions. 
Series of the type (|^),(|) are known to be asymptotic at best 0. We obtain numerical 
values of the exponents Xmn applying Borel transformation improved by a conformal 
mapping procedure p. Tables |l] (H) contain the results for the exponents A^^ i^^rm^) 
describing the scaling of the rate (0) of the reaction (|1|) in the vicinity of a trap attached 
at the core of a polymer star made of m RWs or SAWs respectively. The difference 
between the numbers obtained by resummed e and fixed d -technique may also serve to 
test the accuracy of the data obtained. One can see the decrease of accuracy for high 
functionalities of the star arms caused by the growth of the coefficients of the asymptotic 
expansions (0),(||) with m,n. Let us note, that the case m = 2 corresponds to a trap 
located on the chain polymer, whereas m = 1 corresponds to a trap attached at the 
polymer extremity. 

Table 2: Values of the exponent X^^m^ at d = 3 obtained by e-expansion (e) and by fixed 
dimension technique {3d). 

1^ 1 2 3 4 5 6 

n e 3d e 3d e 3d e 3d e 3d e 3d 



1 


.43 


.45 


.70 


.70 


.89 


.91 


1.03 


1.09 


1.14 


1.24 


1.23 


1.39 


2 


.79 


.81 


1.30 


1.32 


1.69 


1.76 


1.97 


2.12 


2.21 


2.43 


2.41 


2.72 


3 


1.09 


1.09 


1.85 


1.91 


2.41 


2.54 


2.84 


3.10 


3.21 


3.57 


3.53 


4.02 


4 


1.35 


1.37 


2.33 


2.43 


3.07 


3.28 


3.66 


4.02 


4.17 


4.67 


4.60 


5.27 


5 


1.60 


1.64 


2.77 


2.93 


3.69 


3.99 


4.44 


4.92 


5.08 


5.75 


5.63 


6.49 


6 


1.81 


1.89 


3.18 


3.40 


4.26 


4.66 


5.17 


5.77 


5.95 


6.78 


6.62 


7.68 



Numerical values (tables |l|, |) and analytic expressions (||), (j|) for the exponents 
governing the reaction rate in diffusion-controlled reactions involving star polymers are 
the main result of the present study. To conclude, let us analyze several particular cases: 

• For a given m-star absorber oi size R the reaction rate (0) scales as k^n ~ (-R//)"^™". 
Increase of the size i? by a factor of a results in fc^^ ~ {aR/l)~^"^'' leading to: 

k'mn/kmn ~ a"^™ (7) 

with Xmn positive: increasing i? by a factor of a the reaction rate decreases a~'^™" 
times. 

• For a given reaction type (|T|) (i.e. for a fixed number n of particles which are 
trapped simultaneously) attaching mi additional arms to an ?n-arm star absorber 
results in a decrease of the reaction rate: 

kim+rm),n/kmn ~ (i?//)"(^('"+™l)."-^™) , (8) 

as far as A^an > A^in for m2 > mi. 



• For a given m-star absorber the change of the type of reaction (||) to A^^ + B ^ 
0, ni > n results in a decrease of the reaction rate: 

km,in+n,)/k^n ~ (i?//)-(^-.("+"l)-^™), (9) 

as far as Xmn2 > ^mm for n2 > ni. 
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